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The spin structure functions gi for the proton and the deuteron have been measured over a wide 
kinematic range in x and Q 2 using 1.6 and 5.7 GeV longitudinally polarized electrons incident upon 
polarized NH3 and ND3 targets at Jefferson Lab. Scattered electrons were detected in the CEBAF 
Large Acceptance Spectrometer, for 0.05 < Q 2 < 5 GeV 2 and W < 3 GeV. The first moments of g\ 
for the proton and deuteron are presented - both have a negative slope at low Q 2 , as predicted by 
the extended Gerasimov-Drell-Hearn sum rule. The first extraction of the generalized forward spin 
polarizability of the proton 7q is also reported. This quantity shows strong Q 2 dependence at low 
Q 2 . Our analysis of the Q 2 evolution of the first moment of gi shows agreement in leading order 
with Heavy Baryon Chiral Perturbation Theory. However, a significant discrepancy is observed 
between the 7q data and Chiral Perturbation calculations for even at the lowest Q 2 . 
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Fundamental to our understanding of nuclear matter 
is a complete picture of the spin structure of the nu- 
cleon. The spin of the nucleon arises from the spin and 
orbital angular momenta of both the quarks and gluons. 
One way to access the quark spins in lepton scattering 
is through measurements of the spin structure functions 
gi and gi which are not well known at low momen- 
tum transfer to the target nucleon (Q 2 < 2 GeV 2 ). At 
larger momentum transfer, gi(x,Q 2 ) — ^J^efAqi(x) (in 
the parton picture), where Aqijqi is the net helicity of 
quarks of flavor i in the direction of the (longitudinally 
polarized) nucleon spin, qt is the probability of finding a 
quark of flavor i with momentum fraction x, and e, is the 
quark charge. (The Bjorken scaling variable x — -^-^ in 
the lab frame, M is the nucleon mass and v is the energy 
transferred from the electron to the target nucleon.) At 
sufficiently small Q 2 , g\ and its moments can be more 
economically described by hadronic degrees of freedom 
and effective low-energy approximations to QCD, like 
Chiral Perturbation Theory (%PT). 



There is particular interest in the first moment of gi, 
Ti(Q 2 ) = J gi(x, Q 2 )dx, which is related to the frac- 
tion of the nucleon spin carried by quark spins. The 
upper limit of the integral, xq, corresponds to pion pro- 
duction threshold. This limit excludes elastic scattering, 
which otherwise dominates the low Q 2 behavior of the 
integral. Ti is constrained as Q 2 — > by the Gerasimov- 
Drell-Hearn (GDH) sum rule [Id to be -^-Q 2 , where 
k is the anomalous magnetic moment of the nucleon. At 
high Q 2 , Ti has been measured in deep inelastic scatter- 
ing (DIS) experiments at SLACJl i], CERN @, 0, 1 
and DESY [9] . Ji and Osborne [lfj] have shown that the 
GDH sum rule can be generalized to all Q 2 via 

r 1 (Q 2 ) = ^-s 1 (v = o,Q 2 )-rf(Q 2 ), (i) 

where S\(y, Q ) is the spin-dependent virtual photon 
Compton amplitude and T± is the contribution to the 
integral from elastic scattering. At high Q 2 , Si can be 
calculated using the operator product expansion (OPE). 
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By comparing the OPE twist series with rv one can ex- 

" II, El, 0, El, which 
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tract higher twist parameters 
are sensitive to quark-gluon and quark-quark correlations 
in the nucleon at moderate Q 2 . Lattice QCD calculations 
may eventually be available in the moderate Q 2 region 
below the range of applicability of the OPE. At low Q 2 , 
S\ can be calculated in xPT, a model-independent ef- 
fective field theory but it is not clear how high in 
Q 2 these calculations can be applied 17,E3|- Thus Pi 
presents a calculable observable that spans the entire en- 
ergy range from fundamental degrees of freedom (quarks 
and gluons) to effective ones (hadrons). 

Higher moments of g\ are interesting as well. In our 
kinematic domain, these moments emphasize the reso- 
nance region over DIS kinematics because of extra fac- 
tors of x in the integrand. The fundamental generalized 
forward spin polarizability of the nucleon is given by [l9| 



70 (Q 2 ) =C(Q 2 )/ x 2 \ 9l (x, 



UP 

~Q 2 



x 2 g 2 (x 1 Q 2 ) | dx, 

where the kinematic factor C(Q 2 ) — WaM 2 / Q 6 and a is 
the fine structure constant. At high Q 2 one would expect 
cj2 to diminish significantly and g\ to vary logarithmically 
with Q 2 , thus 7q weig hted by Q 6 should be largely inde- 
pendent of Q 2 [3, Il9l . |20T ] . A measurement of 70 on the 
neutron indicates no evidence for such "scaling" below 
Q 2 = 1 GeV 2 , and furthermore the data barely agree 



2l|. No measurement 



with xPT calculations at low Q 2 
of 70 on the proton has been reported so far. 

In order to advance our theoretical understanding of 
the nucleon spin, it is essential to have data on the spin 
structure functions at low Q 2 and in the resonance region, 
as well as at DIS kinematics. Data in the resonance re- 
gion are necessary to calculate moments, especially at low 
and moderate Q 2 . Until recently, data in the resonance 
region were quite scarce 22|. but new measurements of 
spin structure functions in the resonance region have now 
been reported on proton 23|, 24], deuteron 25 1 and 3 He 
targets [26|, [StJ from Jefferson Lab. Testing xPT at low 
Q 2 has increasingly been a focus of new spin structure 
experiments 28l 23. 30j|. 

The EG1 experiment of the CLAS collaboration has 
collected new data using longitudinally polarized 1.6 and 
5.7 GeV electrons on proton (NH 3 ) and deuteron (ND 3 ) 
targets 31] . These data cover a wide kinematic range 
that includes invariant mass W 2 = M 2 + 2Mv — Q 2 from 
elastic scattering (quasielastic for the deuteron) up to 9 
GeV 2 [32I [33| . First results for the generalized forward 
spin polarizability of the proton and new results for the 
first moments of g\ and gf at low and intermediate Q 2 
in the range 0.05 < Q 2 < 3 GeV 2 are reported in this 
letter. 

In the EG1 experiment, the beam was produced from 
a strained GaAs wafer and had an average polarization 
of 70% as measured by Moller polarimetry [34j . The po- 



larization of the electrons was flipped at 30 Hz pseudo- 
randomly. The beam was rastered over the face of the 
target cell to avoid heating and depolarization. The cur- 
rent varied from 0.3 nA to 10 nA depending on the beam 
conditions and target. 

The product of the beam and target polarizations PbPt 
was determined from the data through comparison with 
the known elastic scattering asymmetry and ranged from 
0.50 to 0.60 for the NH 3 target and from 0.12 to 0.23 for 
the ND 3 target. Data were also taken with 12 C, 4 He and 
frozen 15 N to determine the dilution from unpolarized 
materials 35] . 

Scattered electrons were detected in the CEBAF Large 
Acceptance Spectrometer (CLAS) in Hall B [34], covering 
a range in polar angle from 8° to 45° in the efficient 
region of the detector. Data acquisition was triggered 
by a coincidence between the Cerenkov detector and the 
calorimeter in any one of the six sectors. Only electrons 
detected in a region of the Cerenkov detector with an 
efficiency greater than 80% were used in the analysis. 
Additional details about the experiment can be found in 
Refs. ' 



25, 3 



We measured the raw inclusive double spin asymme- 
try with longitudinally polarized beam and target in each 
Q 2 and W bin. This raw asymmetry was then corrected 
for the difference in accumulated charge in the two beam 
polarization states, e + e~ pair production and pion con- 
tamination. Polarization and dilution factors were di- 
vided out and radiative corrections applied. The result- 
ing asymmetry, A/ / , is proportional to a linear combina- 
tion of the two virtual photon asymmetries A\ and A2 
[251 ] . Using a parameterization of the world data to model 
Ai [25] ]. the unpolarized structure function F\ 3(| 37], 
and the ratio of transverse to longitudinal structure func- 
tions R [3^, Ai and g\ were extracted using: 



A^x.Q 2 ) = —A/f-rjAa, 



gi{x,Q 2 



Fi 



l + r 



:{Ai+jA 2 ), 



(3) 
(4) 



where the depolarization factor D depends on R, rj is a 
kinematical factor and 7 2 = . The generalized forward 
spin polarizability for the proton was calculated from the 
data for Ai and the F% parameterization using 70 (Q 2 ) = 
C J* AiFix 2 dx, which is equivalent to Eq. (2). 

The total systematic error on g\ varies greatly depend- 
ing on the kinematic bin; for the proton it is roughly 10% 
and for the deuteron it is typically 15% for the 1.6 GeV 
data and 20% for the high energy data. The systematic 
error is dominated by model uncertainties on Ai , F\ and 
R, which are estimated by using different parametriza- 
tions of the world data. For the deuteron data the un- 
certainty in PbPt also contributes substantially to the 
systematic error. 

The values of g\ and gf were extracted for Q 2 from 
0.05 to 5 GeV 2 and for x greater than 0.1; all results are 
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Q 2 (GeV/c) 2 Q 2 (GeV/c) 2 



FIG. 1: as a function of Q 2 . The data reported here 
(EGlb) are shown as the solid cirlces, along with the ear- 
lier EG1 data (EGla)|U, SLAC [H] and Hermes data [|, 
shown for comparison. The filled circles represent the present 
data, including an extrapolation over the unmeasured part of 
the x spectrum using a model of world data. Phenomeno- 
logical models of Burkert and Ioffe [39l . |40|] and Soffer and 
Teryaev [4lJ are represented by solid and dashed lines, re- 
spectively. The grey band represents the systematic error. In 
the right plot, the scales are expanded and x?T calculations 
from Bernard [131 and Ji 1131 are included. 



available from the CLAS database [38|]. At low Q 2 , the 
A(1232) resonance is quite prominent, with a negative 
asymmetry as expected for this transition. It decreases 
steadily in strength as Q 2 increases. In the mass region 
above the A(1232) resonance, g\ increases from nearly 
zero to large positive values as Q 2 increases. In the 
A(1232) region and at low Q 2 , gf/2 is consistent with 
g\, as expected for a transition to an isospin | state. 
However, at high Q 2 , g^ is significantly larger than gf/2, 
indicating a negative contribution from the neutron. 

The first moments of g\ and gf are shown in Figs. [1] 
and [2l respectively. The parametrization of world data 
[HI is used to include the unmeasured contribution to 
the integral down to x = 0.001. The systematic un- 
certainty (shown by the grey bands) includes the model 
uncertainty from the extrapolation to the unmeasured 
region. Only the Q 2 bins in which the measured part 
(summed absolute value of the integrand) constitutes at 
least 50% of the total integral are shown. For the proton, 
the parametrization is also used at high x (in the range 
1.09 < W < 1.14 (1.15) GeV for the 1.6 (5.7) GeV data). 
For the deuteron, the integration is carried out up to the 



FIG. 2: rf/2 as a function of Q 2 . The symbols and curves 
are the same as for Fig. [T] 



nucleon pion production threshold at high x, excluding 
the quasi-elastic and electro-disintegration contributions. 
Our low Q 2 coverage allows us to observe, for the first 
time, the slope changing sign at low Q 2 , consistent with 
the expectation of a negative slope given by the GDH 
sum rule at very low Q 2 . In general the data are well de- 
scribed by the phenomenological models of Burkert and 
Ioffe \M,\M and Soffer and Teryaev [HI. 

The low Q 2 Ti data are shown in more detail in 
the right-hand panels of Figs. [T] and [2j It is possible 
to make a quantitative comparison between our results 
for and Tf at low Q 2 and the next-to-leading or- 
der xPT calculation by Ji, Kao and Osborne |18|, who 

find rf (Q 2 ) = -^Q 2 + 3.89Q 4 + • • • and T?(Q 2 ) = 

2 

— ■jtjfjQ 2 + 3. 15Q + • ■ • ■ Treating the deuteron as the 
incoherent sum of a proton and a neutron, and correcting 
for the D-state as discussed in Ref. [42j ]. 

r?(Q 2 ) - \{i i.5u; D ) {n(Q 2 ) + r"(<2 2 )} , (5) 

where u>d = 0.056 is the weight of the D-wave in the 
deuteron, one finds that Tf(Q 2 ) = -0.451Q 2 + 3.26Q 4 . 
In the range of Q 2 from to 0.3 GeV 2 , we fit T\ and Tf 
to a function of the form aQ 2 + bQ A + cQ 6 + dQ 8 , where 
a is fixed at —0.456 (proton) and —0.451 (deuteron) by 
the GDH sum rule. Note that the GDH sum rule on the 
deuteron here excludes the two-body breakup part, which 
otherwise nearly cancels the inelastic contribution [43I ]. 
The fit results for the proton, b = 4.31±0.31 (stat) ±1.36 
(syst), and for the deuteron, b = 3. 19 ±0.44 (stat) ±0.68 
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(syst), are both consistent with the Q 4 term predicted 
by Ji et al. The fit (labelled "Poly Fit") is shown 

in the right-hand panels of Figs. [1] and [5] along with Ji's 
prediction. Clearly the Q e term becomes important even 
below Q 2 = 0.1 GeV 2 and this term needs to be included 
in the xPT calculations in order to extend the range of 
their validity beyond roughly Q 2 = 0.06 GeV 2 . The xPT 
4 th order (one-loop), relativistic calculation by Bernard 



et al. 17j is also shown in Figs. Q] and [2] Not shown is the 



result from Bernard et al. that includes an estimate of the 
A(1232) and vector meson degrees of freedom, which are 
important at low Q 2 . That result has large uncertainties 
and is consistent with our data. 



^ 3 

§ 2 
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FIG. 3: Generalized forward spin polarizability 7q as a func- 
tion of Q 2 for the full integral (closed circles), the measured 
portion of the integral (open circles) and Q 2 — [44| (trian- 
gle). The systematic error on the measured (grey) and un- 
measured (dark) contributions are indicated by bands. vPT 
calculations [13,111 are shown along with MAID 2003 Q. 
The data shown on the right are weighted by Q 6 /(16«M 2 ). 
Our parametrization of world data is also shown at moderate 
to high Q 2 . 



Fig.[3]shows the result for the generalized forward spin 
polarizability of the proton 7q (Q 2 )- Since 70 is weighted 
by an additional factor of x 2 compared to T\, the inte- 
gral is mostly saturated by the A(1232) resonance and 
uncertainties due to the low- a; extrapolation are greatly 
reduced. The MAID 2003 model follows the trend of 
the data but lies systematically below them. The MAID 
model is consistent with our data for A\ in the A res- 
onance region, but MAID includes only single-pion pro- 
duction channels, which leads to an underestimation of 
the unpolarized structure function F± entering the defi- 



nition of 70. 

Unlike Ti, 70 is not constrained at Q 2 — and is there- 
fore a more stringent test of Chiral Perturbation calcula- 
tions. The leading order heavy baryon xPT calculation 
by Kao, Spitzenberg and Vanderhaeghen [ZBj], shown by 
the dotted line in Fig. [3j includes the A resonance con- 
tribution. Their 4 th order calculation (dashed line) is 
of opposite sign and shows no sign of convergence; nei- 
ther calculation reproduces the trend or magnitude of the 
data. The relativistic x?T calculation of Bernard, Hem- 
mert and Meissner converges better at 4 order 17 1. 



That calculation, including the resonance contribution, 
is represented by the grey band in Fig. [3l and is also in 
serious disagreement with the data. The A(1232) and 
vector meson contribution is negative (around —2 x 10~ 4 
fm 4 ) and is consistent with the calculation by Kao et al. 
at Q 2 — 0, suggesting that the discrepancy at low Q 2 is 
mainly due to the non-resonance terms [471 ]. 

In the right-hand panel of Fig. [3J 7g is weighted by a 
factor of Q 6 /(16aM 2 ). In the limit of very large Q 2 , this 
expression converges to the third moment of g\ , 02, which 
is expected to scale approximately in the framework of 
OPE. Our data seem to be leveling off above Q 2 = 1.5 
GeV 2 , but do not go high enough in Q 2 to confirm scaling 
behavior. We also show an evaluation of J Q A\F\x 2 dx 
(shaded band) based on our model for F\ and a fit to 
the world data for A\ (mostly from SLAC, SMC, and 
HERMES in addition to our own data). The width of 
the shaded band indicates the combined one-sigma un- 
certainty of the models of F\ and A\. Our model con- 
firms the leveling-off around Q 2 — 2 GeV 2 and shows a 
logarithmic fall-off at higher Q 2 . 

In summary, gi(x, Q 2 ) for the proton and the deuteron 
have been measured over a vastly expanded kinematic 
range at low and intermediate momentum transfer, which 
includes the entire resonance region and part of the DIS 
regime. These measurements enable us to evaluate mo- 
ments of <?i over a wider range in Q 2 , decreasing extrap- 
olation uncertainties. The first extraction of 7^ has been 
reported along with a new precise mapping of and Tf 
down to lower Q 2 than previously available. At mod- 
erately high Q 2 our data for Q 6 7g seem to level off, in 
agreement with models and QCD expectations, and we 
see the expected trend toward DIS results in Ti. It will 
be interesting to extend these measurements to higher 
Q 2 once the upgraded beam energy is available at Jef- 
ferson Lab. At low Q 2 , the first moments of g\ and gf 
exhibit a change in the sign of the slope, to match the 
negative slope constraint from the generalized GDH sum 
rule, and are consistent with xPT calculations for mo- 
mentum transfer values up to about 0.06 GeV 2 . It is 
important to note, however, that these xPT calculations 
also assume the validity of the GDH sum rule; a more 
sensitive test of xPT calculations is 70 (Q 2 )- We observe 
that xPT calculations fail to describe our results for yS, 
even for Q 2 as low as 0.05 GeV 2 . The xPT calculations 
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are increasingly being used to extract results from lattice 
QCD and it is critical to understand their range of ap- 
plicability [lq . Data for the isoscalar quantity 7q — 7q 
have also been published by our collaboration and may 
give additional guidance to future theoretical work in this 
area [IH . We also look forward to results from new exper- 
iments at Jefferson Lab, in which spin structure functions 
down to Q 2 = 0.01 GeV 2 will provide a more stringent 

test of xpt [Uliillaii 
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